DISTURBANCE ESTIMATED TYPE CONTROL SYSTEM, GAS COMPRESSOR CONTROL 
SYSTEM AND METHOD OF DESIGNING A DISTURBANCE ESTIMATED TYPE CONTROL 
SYSTEM 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 
The present indention concerns a disturbance estimated type 

control system, a gds compressor control system and a method of 
designing a disturbance estimated type control system and, in 
particular, relates Jto a disturbance estimated type control system 
and a gas compressor control system that are insensitive to 
parameter variations such as a heat load variation and a rotating 
speed change and a method of designing a disturbance estimated 
type control systen. 

2. Description of the Related Art 

Displacement control of a refrigeration air conditioning 
system for a car air conditioner is for observing an air temperature 
at an outlet of an evaporator or a room temperature by a temperature 
sensor and changing a displacement of a gas compressor such that 
the observed temperature coincides with a target temperature set 
in advance. 

Alternatively, the displacement of the gas compressor is 
changed such that a refrigerant suction pressure of the gas 
compressor becomes a desired pressure value. For example, in 



Japanese Patent No. 1746774, displacement control is performed 
using rotating speed information of a gas compressor together with 
air temperature information. 

An overall schematic configuration view of the air 
conditioning system is shown in Fig. 12. The air conditioning 
system is mounted, for example, on a car. In Fig. 12, an evaporator 
51 is a heat exchanger involving ventilation by a fan 52 and cools 
air in a cabin. Refrigerant gas that is evaporated within the 
H evaporator 51 is taken into a variable displacement type gas 

s * 

a compressor 10 and pressurized to be sent to a condenser 53. 

fu . 

Q In the condenser 53, the refrigerant gas is liquidized and 

^ heat absorbed from the inside of the cabin is discharged to the 

IS 

\± outside of the car. An expansion valve 54 is made to rapidly 

ru 

nj decrease a pressure of the refrigerant gas from a high pressure 
□ to a low pressure. A rotating shaft 11 of the variable displacement 
type gas compressor 10 is made such that a shaft power of an engine 
59 is transmitted to it and it is rotationally driven. 

A sectional view of the variable displacement type gas 
compressor 10 is shown in Fig. 13 and a sectional view cut along 
line A-A and viewed in the direction of arrows A in Fig. 13 of 
the variable displacement type gas compressor 10 is shown in Fig. 
14. 



A suction port 1 of the variable displacement type gas 
compressor io is made to suction refrigerant gas from the evaporator 
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51 connected to the outside. 

A cylinder 3 is clamped between a front head 5 and a rear 
side block 7. A rotor 9 is rotatably arranged within the cylinder 
3. 

The rotor 9 is penetrated by and fixed to a rotary shaft 
11. Bane grooves 13 are formed in the radial direction on the 
external circumference of the rotor 9 and banes 15 are slidably 
inserted in the bane grooves 13. Then, the banes 15 are biased 
to the internal wall of the cylinder 3 by a centrifugal force and 
an oil pressure of the bottoms of the bane grooves 13 when the 



ru 

/■* rotor 9 is rotating 



*y The inside of the cylinder 3 is partitioned into a plurality 



of small chambers by the rotor 9 and the banes 15. These small 

U 

y chambers are referred to as compressing chambers 17 , which repeat 
□ increasing and decreasing changes of a displacement by the rotation 
of the rotor 9. 

Then, when the rotor 9 rotates to change the displacement 
of the compressing chambers 17 in this way, low-pressure 
refrigerant gas is taken in from the suction port 1 and compressed 
by the displacement change. A case 19 is fixed in the peripheral 
end portion of the cylinder 3 and the rear side block 7 and a discharge 
chamber 21 is formed inside this case 19. 

High-pressure refrigerant gas compressed in the compressing 
chamber 17 is sent to the discharge chamber 21 via a discharge 
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port 23 and a discharge valve 25. Then, the refrigerant gas is 
sent to the condenser 53 in the outside from the discharge chamber 
21 via the discharge port 27. 

This variable displacement type gas compressor 10 is 
provided with a displacement varying mechanism 30. This 
displacement varying mechanism 30 is made such that it can variably 
adjust a discharge capacity of refrigerant gas according to a 
temperature in a cabin. An example of a configuration of the 

N* displacement varying mechanism 3 0 is shown in Fiq. 15. 

□ 

P } A control plate 29 is arranged in the front head 5 to face 

y the side portion of the cylinder 3. Notches 29a are provided in 

CP 

yy two parts of the control plate 29. These notches 2 9a cause the 

3 

m inside of the cylinder 3 and a suction chamber 31 leading to the 

fU 

fU suction port 1 to communicate with each other. On the other hand, 

m 

Q the compressing chambers 17 are formed in spaces closed by the 
part on the control plate 2 9 where there is no notch, the internal 
wall of the cylinder 3 and the banes 15. 

^l^^ When the control plal e 2 9 is rotated to the right , the notches 

direction, whereby the positions where 
are formed also move to the right side 
compressing chambers 17 at this point 
is also reduced. In this way, a discharge capacity is adjustable 
by rotating the control plate 29. 

The rotation of the control plate 2 9 is performed by a driving 



29a is rotated in the right 
the compressing chambers 1 
and the displacement of the 
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shaft 39 of hydraulic drive via a pin 33. Oil is injected into 
a sleeve 35 from the discharge chamber 21 by adjusting the opening 
degree of a control valve 3 7 and the driving shaft 3 9 is moved 
straight by oil pressure at this point. Then, this straight 
movement is converted into rotating movement via the pin 33 to 
cause the control plate 29 to rotate. 

An injection amount of oil can be changed by altering the 
opening degree of the control valve 37. This alteration of the 
opening degree is performed by changing a displacement control 
command value (duty ratio) shown in Fig. 16. 

The control plate 2 9 is rotated under the balance with an 
elastic force by a spring 38 in accordance with a differential 
pressure between a control pressurePc in the sleeve 35 and a pressure 
Ps in the suction chamber 31. 

Further, in Fig. 12, for example, in order to detect an 
air temperature at the outlet of the evaporator 51, a temperature 
sensor 55 is disposed. In addition, in order to detect a rotating 
speed of the engine 59, a rotating speed sensor 57 is disposed. 

Then, a displacement control command value is calculated 
in a control circuit 61 based on a detected signal of this temperature 
sensor 55 and a rotating speed of the engine 59 . This displacement 
control command value is signal-amplified by a displacement control 
signal generating circuit 65 and, then, transmitted to the control 
valve 37 of the displacement varying mechanism 30. 
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Next, a control method of this displacement varying 
mechanism 30 will be described based on a flow chart of Fig. 17. 

For simplicity, description will be made with the case in 
which a displacement of the variable displacement type gas 
compressor 10 is reduced such as the case in which a rotating speed 
of the engine 59 is increased as an example. 

Now, the case in which a detected temperature at the outlet 
of the evaporator 51 has become lower than a target temperature 
67 at the outlet of the evaporator 51 is assumed. In this case, 
in order to prevent the inside of the cabin from being cooled 
excessively, it is necessary to lower cooling capability. 

First, in step 1 ( it is abbreviated as SI in the figure 
and the same is true for subsequent steps), a target refrigerant 
flow rate of the variable displacement type gas compressor 10 is 
calculated based on a temperature deviation between a target 
temperature and a detected temperature. The calculation in this 
case is performed by PID control or the like. Next, in step 3, 
a rotating speed of the variable displacement type gas compressor 
10 is calculated from this calculated target refrigerant flow rate 
or a discharge capacity of the variable displacement type gas 
compressor 10 is calculated taking a rotating speed of the engine 
59 into account. 

In step 5, a rotating speed correcting calculation is 
performed from this discharge capacity based on a not-shown 



property curve showing a relation between a discharge capacity 
and a displacement control command value and a displacement control 
command value by which an opening degree of the control valve 3 7 
should be adjusted is determined. The displacement control 
command value at this point is commanded to be small . As a result, 
an average current falls in step 7 and the opening degree of the 
control valve 3 7 is reduced in step 9. 

At this point, the control pressure Pc inside the sleeve 
M= 35 falls in step 11. Thus, the driving shaft 39 is moved downward 

2 

p in step 13 and the control plate 29 rotates to the right in step 

. y 15. As a result, the discharge capacity of the variable 

SI 

ry displacement type gas compressor 10 becomes small in step 17 and 

s 

I* the cooling capability falls. 

ru 

?U In addition , other than the above-described control method , 

• m 

O methods such as a method of preparing a plurality of control modes 
and switching among the modes and a method of altering a model 
based on observation information and changing a parameter of a 
compensator according the model alteration have been used because 
the air conditioning system has different physical property 
depending on a heat load variation and a difference of an operation 
condition and also has strong nonlinearity (Japanese Patent Nos. 
2090676 and 2086982) . 

However, the embodiment according to Japanese Patent No. 
1746774 increases the number of sensors because a rotating speed 
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sensor is required. In addition, even if an existing rotating 
speed sensor is used, there is a disadvantage in that the number 
of communication means increases. Thus, although a technology 
for realizing equivalent or higher control performance even without 
rotating speed information is necessary, the realization has been 
difficult conventionally. 

On the other hand, in the case of Japanese Patents Nos. 
2090676 and 2086982 , providing a plurality of control compensators 
M> leads to increase of development load. In addition, alteration 
Q of a compensator is likely to cause trouble in practical use because 

ru 

ri there arise problems in that the number of sensors is increased, 

Ul 

ffi 

p\ development load is increased and the compensator itself becomes 

CO 

rU-x^vO y The present invention has been devised in view of such 
= ^/Conventional problems and iJt is an object of the present invention 
to provide a disturbance estimated type control system and a gas 
compressor control system that are insensitive to parameter 
variations such as a heat load variation and a rotating speed change 
and a method of designing a disturbance estimated type control 
system. 

SUMMARY OF THE INVENTION 
Therefore, according to the present invention, there is 
provided a disturbance estimated type control system, comprising: 
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a control object; disturbance estimating means for estimating a 
disturbance added to an input of the control object based on an 
input signal to be inputted in the control object and a detection 
signal detected from the control object; compensating means for 
compensating for the detection signal and outputting a control 
signal; and calculating means for subtracting a disturbance 
estimated by the disturbance estimating means from the control 
signal and considering a subtraction result as the input signal , 
characterized in that a mathematical model (nominal model) of the 
control object comprising of a state equation and a transfer 
function is prepared by system identification based on an 
experiment , the compensating means is designed from the 
mathematical model ( nominal model ) , and the disturbance estimating 
means is prepared from an expanded state equation (expansion 
system) comprising of a mathematical model ( expansion model ) taking 
a disturbance into account with respect to the mathematical model 
(nominal model). 

A nominal model of a syitem is prepared utilizing a system 
identification technology. Then, disturbance estimating means 
is prepared from am expanded state equation (expansion system) 
comprising of mathematical model (expansion model) taking a 
disturbance into account wit/h respect to this mathematical model 
(nominal model). J 

A disturbance against a control object can be offset by 



ru 



m 
m 



ru 



N 1 



subtracting a disturbance estimated by the disturbance estimating 
means from a control signal that is an output of compensating means . 
Thus, a system control system insensitive to a parameter variation 
caused by a change in physical property can be realized with a 
fixed compensator having a simple structure. 

Further , the present invention is configured by comprising : 
a variable displacement type gas compressor having displacement 
altering means that is capable of altering a displacement within 
a compressing chamber; an input signal to be inputted in the 
displacement altering means; a detection signal in which at least 



SI one piece of information among an air temperature in a room, an 



air temperature at an outlet of an evaporator, a refrigerant flow, 
a refrigerant pressure on a suction side of a gas compressor and 
the like is detected; disturbance estimating means for estimating 
a disturbance of the variable displacement type gas compressor 
based on the detection signal and the input signal; compensating 
means for compensating for the detection signal to output a control 
signal; and calculating means for subtracting a disturbance 
estimated by the disturbance estimating means from the control 
signal to consider a subtraction result the input signal. 

A disturbance is, for example, a variation of a fan 
ventilating an evaporator, a variation of the number of rotations 
of a variable displacement type gas compressor that is driven by 
an engine, or the like. A refrigeration air conditioning system 
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control system, which is insensitive to a parameter variation 
caused by the change in physical property values of air and 
refrigerant due to a heat load variation and the change in 
refrigerant flow due to a variation of the number of rotations , 
can be realized with fixed compensating means having a simple 
structure. 

Further, according to the present invention, there is 
provided with a method of designing a disturbance estimated type 

y, control system that comprises: a control object; disturbance 

Q 

Q estimating means for estimating a disturbance added to an input 

OJ 

of the control object based on an input signal to be inputted in 
the control object and a detection signal detected from the control 



ru 

object; compensating means for compensating for the detection 

ru . 

jy signal and outputting a control signal; and calculating means for 
p subtracting a disturbance estimated by the disturbance estimating 
means from the control signal and considering a subtraction result 
as the input signal, characterized by comprising: preparing a 
mathematical model ( nominal model ) of the control object comprising 
of a state equation and a transfer function by system 
identification; preparing an expanded state equation (expansion 
system) comprising of the mathematical model and a mathematical 
model of a disturbance; designing the disturbance estimating means 
from the state equation (expansion system); and designing the 
compensating means from the mathematical model (nominal model). 
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Further, the present invention is characterized by 
comprising: determining whether or not the state equation 
(expansion system) is observable; and if it is determined to be 
unobservable, compulsorily adding an error of 10% or less to a 
coefficient corresponding to an A matrix and/or a C matrix of the 
state equation (expansion system) or a zero-th dimension term of 
a transfer function numerator and preparing an expansion system 
including the error, thereby establishing observability. 

H i In this way, the disturbance estimating means and the 

Q 

P compensating means can be designed even if the state equation 

fU 

~~i (expansion system) is unobservable. 

m 
m 

L, BRIEF DESCRIPTION OF THE DRAWINGS 

ru 

fU Fig. 1 is a block diagram of an expansion system including 
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a nominal model and a disturbance. 

Fig. 2 is a block diagram of a control object model and 
an observer. 

Fig. 3 is a block diagram of a control system. 

Fig. 4 is flow chart of development. 

Fig. 5 is flow chart of procedure of system identification. 
Fig. 6 is an example of an M sequence signal. 
Fig. 7 is an example of an M sequence data used in 
refrigeration system control. 

Fig. 8 is a relation between a rotating speed with respect 
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to environmental conditions and a displacement ratio at a balanced 
point . 

Fig. 9 is an example of input and output data. 

Fig. 10 is an example of detrended input and output data. 

Fig. 11 is an example of an evaluation result. 

Fig. 12 is a schematic view showing an entire configuration 
of an air conditioning system. 

Fig. 13 is a sectional view of a variable displacement type 
gas compressor. 

Fig. 14 is a sectional view cut along line A-A and viewed 
in the direction of arrows A in Fig. 13. 

Fig. 15 is an example of a configuration of a displacement 
varying mechanism. 

Fig. 16 is a displacement control command values (Duty 
ratios ) . 

Fig. 17 is flow chart of control method of displacement 
varying mechanism. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
Embodiments of the present invention will be hereinafter 
described. 

The present invention is for applying a disturbance observer 
whose effectiveness is confirmed in a mechanical system such as 
a robot and a machine tool (see Kohei Ohnishi, Kiyoshi Ohishi, 
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Kunio Miyaji: "A control method of a separate excitation DC machine 
using a state observer", the Institute of Electrical Engineers 
of Japan Thesis Magazine B104-6 , pp. 373-37 9 , 1984) to displacement 
control of a refrigeration air conditioning system. 

Its characteristic is that it is insensitive to a variation 
despite the fact that it is a compensator of a fixed gain and, 
moreover, a simple compensator of low dimension. However, even 
if it is attempted to apply a technology of another field of 
application to the refrigeration air conditioning system and divert 
it to the field, even a nominal model that becomes a premise of 



J designing a disturbance observer cannot be easily prepared because 
the refrigeration air conditioning system is complicated and has 



jV strong nonlinearity. That is, the disturbance observer technology 

hi 

pj that has been used for a motor cannot be diverted without change. 



Q Thus, the system identification technology is utilized to 

prepare a nominal model of a refrigeration air conditioning system, 
an expansion system comprising of the nominal model and a newly 
assumed disturbance is configured and, after confirming conditions 
for the expansion system to meet observability that is a 
configuration condition of an observer, a new disturbance observer 
is designed to make it possible to mount it on an actual machine. 

Incidentally, the term "nominal" , which means " in name only " , 
is based on the fact that, since an actual system varies, a 
mathematical model is not limited to its mathematical model . Here, 
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it refers to a model that is found according to one average condition . 

First, a method of obtaining a nominal model will be 
described • 

It is assumed that a refrigeration air conditioning system 
is modeled as in the following expressions 1 and 2 by means such 
as system identification. 

[Expression 1] 

[Expression 2 ] 

y = C P n* pn 

However, it is assumed to be an actual vector or an actual 
matrix of x pn (n pn X l ) , u( 1 X l ) ,y ( l X n pn ) , A pn ( n pn X n pn ) , B pn ( n pn X 



M, 1 ) , C pn ( 1 X n pn ) . 

ru 

pj Here, the input u is a driving input signal of an actuator for 



causing a displacement of a variable displacement type gas 
compressor to change and is a voltage in most cases, but may be 
variable data or the like that plays the same role inside a computer 
system for control. 

In addition, the output y is a detected value of an air 
temperature of a part that is to be controlled and is observed 
in the form of a voltage in most cases , but may be variable data 
corresponding to the temperature captured into a computer for 
control . 

Whichever it is, the above-described model is a model under 
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a certain condition and does not include a physical parameter 
variation or a disturbance . This is referred to as a nominal model . 

Next, a method of the system identification will be 
described. 

The system identification per se is a conventional modeling 
theory and technology (see, for example, Shuichi Adachi: System 
identification for control by MATLAB, Tokyo Denki University Press , 
1996). There are known a prediction error method, a maximum 
y> likelihood estimation method, a least-square method, a subspace 

b 

O identification method and the like, and it is assumed that any 

Nl method may be used but representation is finally made in the form 

W 

pi of a state equation. 

In addition, since an identified model is a model that should 

ru 

also have been found by an experiment using input and output, it 

CD 

□ is controllable and observable (see, for example, Kosato, Mita: 
Guide to system control theory, Jikkyo Shuppan Co., Ltd., 1979). 
Conversely, even if there is an uncontrollable or unobservable 
part in an actual system, modeling is performed with that part 
neglected. 

In the case of a one-output system, a necessary and sufficient 
condition of controllability is that a rank of a controllable matrix 

U Cpn ( n pn X n pn ) of the following expression 3 is n pn . 

[Expression 3] 
U Cpn = \B pn ,A pn B pn ,Al n B pn ,-,A n p :- x B p \ 
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Then, rank (Ucpn)=n P n f that is , 



[Expression 4] 



det(£/ Cp „) * 0 



In the case of a one-output system, a necessary and sufficient 
condition of observability is that a rank of an observable matrix 
U 0pn ( n pn X n pn ) of the following expression 5 is n pn . 

[Expression 5] 



observability taking a disturbance into account will be described . 

Correction taking a disturbance into account is added to 
the above-described nominal model. 

A parameter variation and an unknown external input are 
collectively considered a disturbance d at an input end to define 
a new state amount x c i=d and prepare expansion systems of a state 
equation Expressions 7 and 8. 

[Expression 7] 
x = Ax + Bu 



U 0p n - [ C pn> C pn A pn> C pn A ln' " > 



C A""'*] 7 

pn J 



< Then, rank(U 0pn )=n pn , that is, 



[Expression 6 ] 
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[Expression 8] 

y = Cx 

Here, Expression 9 is assumed. 
[Expression 9] 



Here, each element is 
[Expression 10] 















X = 








,B = 










_ 0 . 



A dimension of the system is n=n pn +l at A((n pn +l)X 
(n pn +l) ,B( (n pn +l) XI) ,0(1 X (n p „+l) . 

A block diagram of an expansion system including a nominal 
model and a disturbance is shown in Fig. 1. It is known that, 
if Expressions 7 and 8 are observable, an observer can be designed. 
A necessary and sufficient condition of observability is that, 
in the case of a one-output system, a rank of an observable matrix 

U 0 (nXn) of the following Expression 11 is n. 
[Expression 11] 

U 0 = [C,CA t CA\~,CA"- ] ] r 

That is, Expression 12 must be established. 

[Expression 12 ] 
det(^) * 0 

In the case of this system, an observable matrix is found 
to be Expression 14 from Expression 13. 
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[Expression 13] 
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A B 

pn pn 
Ix n^ 



OA = [C fn 0] 
CA^[C pn A pn C„B„] 



CA 3 = [C pn Al C pn A pn B pn ] 



= \C A C B 1 

^ p« pn pn pn J 

= fcX. C pn A pn B p \ 
= [C pn K C pa Al„B pn ] 



A B 

pn pn 



CA n l = \C A"' 3 C A n A B ] 

L Pft pn pn pn pn J 

CA" ' = \C A n i C A n i B 1 

I pn pn pn pn pn J 

[Expression 14] 



A pn B pn 



^xn„ 

A „n B pn 



= Mr o pn A- 3 B pn ]=[c n A;r c pn A;r'B pn ] 
= [c,„a;: c pn A;; l B pn ]=[c pn A"; n - c pii a"-'b,\ 



U„ = 





0 








C A 

pn pn 


C B 

pn pn 






Y 


C A""'* 

pn pn 


C A n "- 2 B 

pn pn pn 




C A"" 

pn pn 


C A n "' x B 

pn pn pn _ 


C A"" 

pn pn 


C A npa ~*B 

pn pn pn _ 









where Y=\o,C B , -,C A n "' 2 B Y 

[ pn pn r 9 pn pn pn J 



Here, if Expression 15 that is a formula (see Kosato, Mita: 
Guide to system control theory, Jikkyo Shuppan Co., Ltd., 1979) 
is utilized, the determinant of Expression 12 can be represented 
as Expression 16 from Expression 4. 

[Expression 15] 
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det 



X Y 
Z W 



= det(*) • det^ - ZXY) for det(x) * 0 



[Expression 16] 
det(U 0 ) = det(U 0an )det(c A n X B - C A"" • £/"' y) 

\ O J \ Opn / y pn pn pn pn pn Opn J 

If these values are not zero, Expressions 15 and 16 are 
observable, and thus an observer can be designed. 

Next, a pole of an expansion system and controllability 
will be described. 

A pole increased by preparing an expansion system will be 
confirmed. A formula for a determinant (see Kosato, Mita: Guide 
to system control theory, Jikkyo Shuppan Co., Ltd., 1979) 

[Expression 17] 



det 



X Y 
0 W 



= det(^)det(M/) 



is used, from Expression 7 , since the pole of the expression system 
is s satisfying Expression 18, 
[Expression 18] 



det(s/->4)= det 



S hn' A pn -Bpn 



it is seen that it consists of a pole and the origin of an original 
nominal model. 

In the case of a one- input and one-output system, an unchanged 
zero point in the modern control theory and a zero point ( = transfer 
zero point) in the classic control theory coincide with each other 
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and uncontrollability and unobservability occur with a pole and 
a zero point coinciding with each other (see Kosato, Mita: Guide 
to system control theory, Jikkyo Shuppan Co., Ltd., 1979). 

Incidentally, it can be seen without checking a controllable 
matrix that this expansion system is uncontrollable. This is 
because a disturbance is deliberately treated as one of state 
amounts despite the fact that it is physically one of inputs and 
its origin pole cannot be altered by control input as is seen from 
the block diagram in the upper part of Fig. 1. 

It is confirmed that a necessary and sufficient condition 
of controllability is that, in the case of a one-output system, 
a rank of a controllable matrix Uc(nXn) of Expression 19 is n. 



m [Expression 19] 



That is, Expression 20 must be established. 
[Expression 20] 
det(*/ c )*0 

When Expression 2 0 is expanded, Expression 21 is found 
utilizing Expression 17. 
[Expression 21] 



det(U c ) = det 



B A B ••• A""' X B A""B 
0 0 0 0 0 



= det 



^ cpn ^pn ^ pn 



= det(^).0=0 



Therefore, Expression 20 is uncontrollable. In Expression 
18, since the pole of the expansion system consists of a pole and 
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the origin of a nominal system and, moreover, the nominal system 
is controllable, a pole that causes uncontrollability , that is, 
an uncontrollable pole is only the origin 0. 

To supplement the above, since the expansion system is used 
as a tool for designing an observer and is not used for designing 
a compensator, these expansion systems Expressions 7 and 8 may 
be uncontrollable. 

Next, conditions of observability will be described. 
H 5 Since there is a case in which an uncontrollable pole 

s — t 

TBS* 

P simultaneously becomes an unobservable pole, evaluation of 

=--_? 

jj observability is necessary. 
f\ (Case 1) 

s 

M, It is assumed that there is a positive integer q and 

fU • 
pj [Expression 22] 

81 

□ C B = C A B = -=C A"- 2 B =0,C A'-'B *0,q<n 

p pn*^ pn pn *pn*^ pn pn pn pn pn pn pn ' ^ pn 

is established. It is known that, in a one-input and one-output 
system, r=n pn -q represents the number of zero points (Furuta, Kawaji, 
Mita, Hara: Mechanical system control, Ohmsha Ltd., 1984). 

If r=n pn -q=0 by chance in a nominal mode, since there is 
no zero point in the nominal model itself and the Y matrix becomes 
the 0 matrix, the following formula Expression 23 can be used for 
an expansion system (see Kosato, Mita; Guide to system control 
theory, Jikkyo Shuppan Co., Ltd., 1979), 

[Expression 23] 
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= det(*)det(w) 



From Expressions 6, 22 and 23, 
[Expression 24] 
det(U Q ) = det^^ detf^;:"^) * <> 

and the expansion system Expressions 7 and 8 becomes observable. 
(Case 2) 

If r=N pn -q>0, that is, there are one or more zero points 
in a nominal model itself , it is necessary to confirm that an origin 
pole is not turned into an unobservable pole and there is no other 
way but to directly check Expression 16. 

Next, a method of designing a disturbance observer will 
be described. 

An observer is designed using the expansion systems 
Expressions 7 and 8. 

(Design of an identical dimension disturbance observer) 
[Expression 25] 

x = Ax + Bu - K(Cx - y) = (A - KC) x + Bu + Ky 



where x = 





X 




A B 




'B pn ' 




pn 


,A = 


pn pn 












0 . 



x e R n *\A e R nx \B e R nx \C e R U \K eR nx \k pn e R n ^\n = n pn + 1 

and K is selected such that real parts of all peculiar values of 
the A-KC matrix of Expression 25 become negative. A block diagram 
of a control object model and an observer is shown in Fig. 2. In 
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Fig. 2, an input u is, for example, a displacement control command 
value and an output y is a detected temperature of the outlet of 
the evaporator 51. 

(Design of a minimum dimension observer) 

A designing method called a nonsingular form of Gopinath 
may be used for designing the observer. 
(Observer of a digital format) 

The expansion system Expressions 7 and 8 may be discretized 
with a zero-th dimension hold, as a premise and designed as an 
identical dimension observer of a digital format or a minimum 
dimension observer. 

Next, utilization in a control system will be described. 

A block diagram of a control system is shown in Fig. 3. 
A disturbance 103 is applied to a control object 101. An input 
signal 105 (e.g., a displacement control command value) and an 
output signal 107 (e.g., a detected temperature at the outlet of 
the evaporator 51 ) are inputted in an observer 109 and a disturbance 
is estimated. 

A deviation between the output signal 107 and a target value 
111 is calculated by a subtractor 113. This deviation is subject 
to compensations such as a state feedback, a servo system of a 
1 type and a PID control by a compensator 115 and outputted as 
a control signal. This control signal is subtracted from a 
disturbance estimated value from the observer 109 by a subtractor 
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117. 

Therefore, the disturbance estimated value among estimated 
values of the observer is fed back to cancel an actual disturbance. 
In addition, compensations such as a state feedback, a servo system 
of a 1 type and a PID control are performed using another state 
estimated value. 

Next, measures at the unobservable time will be described. 

First, if the above descriptions are arranged as a 
development procedure, it is as shown in Fig. 4. In a development 

s 

O flow of Fig . 4 , a nominal model of a refrigeration air conditioning 

ru 

^ system is prepared utilizing a system identification technology 
in step 31. In step 33, an expansion system is configured from 
^ the nominal model and a disturbance that is assumed anew. 

ru 

ry Next, in step 35, it is determined whether or not the 

p expansion system satisfies observability that is a configuration 
condition of the observer. If it satisfies observability, a new 
disturbance observer is designed in step 37. This disturbance 
observer and a compensator are combined in step 39 and mounted 
on an actual machine in step 41. 

In step 35, if the expansion system does not satisfy 
observability, unobservability is avoided by deliberately adding 
a deviation to an identification model. As a method, a constant 
term of a small value is added to deviate a zero point. 

In this way, a refrigeration air conditioning systemcontrol 
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system, which is insensitive to a parameter variation caused by 
a change in a physical property value of air and refrigerant due 
to a heat load variation or a change in a refrigerant flow due 
to a change in a rotating speed, can be realized with a fixed 
compensator having a simple structure. 

Further, a specific method of system identification will 
be described for reference. 

System identification means preparing some "mathematical 

M> model" that can prove that a system is "identical" with an object 

□ 

^ under a certain "purpose" from a measured value of input and output 

^ data of a dynamic system to be an object. System identification 

01 

pj intended to be performed in the present invention is "system 

s 

1^ identification for the purpose of control" and several methods 

ru 

jfjj have been established. 

01 

O A procedure for system identification follows a flow of 

Fig. 5. In step 51, M sequence data is prepared. Next, in step 
53, implementation in a controller is performed and, in step 55, 
measurement is performed and the data is stored. In step 57, an 
identification program is prepared and, in step 59, identification 
is executed. More specifically, A pn , B pn and C pn are found. 

Then, in step 61, the found A pn , B pn and C pn are evaluated. 
If an evaluation result is displeasing, the identification program 
is revised and the identification work is executed again to perform 
evaluation again. In addition, it becomes necessary to repeatedly 
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perform identification with respect to a control object with 
significantly changing environmental conditions by changing an 
amplitude, a cycle and a length of the M sequence data in various 
ways . Next , preparation of the M sequence data will be described • 
For system identification, Pseudo Random Binary Signals 
(PRBSs) are inputted in a control object and an output is observed. 
Among the signals to be inputted , most frequently utilized signal 
is an M sequence signal (Maximum-length linear shift register 

M, sequence). An example is shown in Fig. 6. 

Q 

Q Next, a specification of the M sequence signal in the case 

rU . 

^ in which it is actually used for refrigeration system control will 

W 

be considered. An example of the M sequence data used for 
refrigeration system control is shown in Fig. 7. Fig. 7 shows 

ru 

pj data strings in which 0 to 10 are command voltages corresponding 

2 i "» 

p to a driving current zero to maximum current of the displacement 

a . 

varying mechanism 30 and is an M sequence signal deflected upward 
and downward by an amplitude 2 with a heat load as a medium load 
and with an average command value 7.85 at the time when a rotating 
speed of the variable displacement type gas compressor 10 is 800 
rpm as an average value. 

As the average value, a balanced point of inputs at the 
time when a target temperature coincides with an air temperature 
under a condition that a heat load and a rotating speed of the 
variable displacement type gas compressor 10 are regulated is used . 
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Since the amplitude depends on nonlinear evaluation, increasing 
and decreasing of the amplitude are attempted in a try and error 
manner. However, the amplitude of less than zero or more than 
ten is unacceptable because it causes saturation, the average value 
is deviated so as to avoid saturation. 

Fig. 8 shows a relation between a displacement ratio (= 
Duty 2) and a rotating speed at the balanced point in five kinds 
of environmental conditions . 

§=* In a temperature control system, this value may be used 

Q 

O as a feed forward value corresponding to a rotating speed. 

ru 

M Incidentally, 0 to 10 used in the above-described M sequence 

UJ ' 



fU 



ru 



correspond to the 0 to 1 . Thus, for example, since a value at 
the time of operation at a rotating speed 800 rpm under the 
environmental conditions of (1) is approximately 0.7, it is 

m 

p sufficient to use ,7 that is ten times as large as the value as 
the average value. Besides, if temperature control is attained 
for the time being, a balanced point at the time of control may 
be used. 

Next, implementation in a controller and execution will 
be described. 

The controller outputs an M sequence signal as a voltage. 
Then, this M sequence signal is given as a control input to the 
displacement varying mechanism 30 of the variable displacement 
type gas compressor 10 that is a control object. A resulting air 
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ru 
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CP 
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m 




temperature at an outlet of an evaporator is detected in the form 
of a voltage and captured in the controller. 

Results of input and output in an identification experiment 
using a control bench are shown in Fig. 9. On the other hand, 
an average value of input and output is eliminated as shown in 
Fig. 10. That is, an average is zero. System identification 
calculation is performed based on this result and A, B and C matrixes 
of a state equation are obtained. 

Next, as shown in Fig. 11 as an evaluation result, validity 
of an obtained mathematical model can be confirmed by comparing 
an actual output and an output at the time when an identical M 
sequence signal input is given to an identical model. 

As reference documents, there are the following. 

1) L. Ljung: System Identification Tool box User's Guide, Math 
Works (translation supervised by Shuichi Adachi) 

2) Shuichi Adachi: System identification for control by MATLAB, 
Tokyo Denki University Press, 1996 

3) Kiyoshi Wada: Commentary "What is the partial space 
identification method?" , Keisoku to Seigyo 36-8 , PP. 569-674, 1997 

As described above, according to the present invention, 
r a system control system that is insensitive to a parameter variation 
caused by a change of a physical! property value or the like can 
be realized with a fixed compensator having a simple structure. 
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